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We study the question of generalizing light-front field theories to finite temperature. We show that 
the naive generalization has serious problems and we identify the source of the difficulty. We provide 
a proper generalization of these theories to finite temperature based on a relativistic description of 
thermal field theories, both in the real and the imaginary time formalisms. Various issues associated 
with scalar and fermion theories, such as non-analyticity of self-energy, tensor decomposition are 
discussed in detail. 
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I. INTRODUCTION 



Quantum field theories are conventionally quantized on a space-like surface. However, quantization on a light-like 
^ . surface provides an interesting alternative which dates back to the works of Dirac . Light-front field theories ^ , 
namely theories quantized on a light-front, have found applications in various branches of physics such as QCD, string 
, theories and membrane theories, among others Q]. In the case of QCD, for example, they provide a method for 
studying non-perturbative phenomena systematically, which is different from the usual lattice studies. 
^-H I In fact, there is a lot of activity in the study of light-front field theories within the context of QCD. Many of 
^ the features of light-front theories are quite distinct from the conventional equal-time field theories - one of the most 
\ significant being that light-front theories are first order (as opposed to conventional theories which are second order) in 
I time derivatives and, correspondingly, describe different degrees of freedom. While the behavior of conventional field 
theories at finite temperature are quite well understood by now 1^, ^, 0] , a systematic study of the thermal properties 
of light- front field theories is lacking so far. It is the purpose of this paper to work out the essential properties of 
' scalar and fermion light-front field theories at finite temperature. We will defer the study of light-front gauge theories 
, to a future publication. 

■ It would seem that the generalization of the thermal field concepts to light-front theories should be straightforward. 
In fact, in n-dimensions, if we define the light-front time variable as 

Oh' 1 

D, x+ ^ ^ (x° + x''-') (1) 



then, the Hamiltonian of the system can be identified with 



^ : H = p-^^{p"- (2) 

This would suggest that the ensemble average of an arbitrary operator, O, in the light-front field theory, should be 
defined as 

{0)[3 = Tr e-l^" = Tv e^^^^ O (3) 

where (3 represents the inverse temperature in units where the Boltzmann constant is unity. This has been, in fact, 
the general thinking However, as we will show, this straightforward generalization is incorrect and leads to various 
problems. We will discuss the proper description of thermal field theories (for scalars and fermions), on the light-front, 
both in the imaginary as well as the real time formalisms. In section II, we will describe the naive generalization 
of the concepts of thermal field theories to light-front scalar field theories, both in the imaginary and the real time 
formalisms, and show that this leads to various problems. The source of the difficulty is identified in section III, where 
we give the proper description of light-front scalar theories at finite temperature. We calculate various quantities of 
interest such as the thermal mass correction as well as the non-analyticity in the self-energy. We show that there 
are more possible limits that can arise in these theories, in contrast to conventional thermal field theories. A short 
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description of the tensor decomposition, which has a richer structure in such theories, is also given. In section IV, 
we discuss briefly the generalization of light-front fermion theories to finite temperature. Here a new feature arises 
since these theories have only half the number of independent degrees of freedom |^ . The generalization of light-front 
gauge theories to finite temperature as well as various other applications are under study and will be described in a 
future publication. 

II. NAIVE GENERALIZATION TO FINITE TEMPERATURE 

In this section, we will discuss the naive generalization of the techniques of thermal field theory to the light-front 
scalar field theories. Let us briefly establish the notation. In n-dimensions, we define 

x± = i=(a;"±x"-i) (4) 

where x"*" is identified with the light-front time coordinate. Denoting the coordinate vector as = (a;"'', x~ ^x), where 
X represent the {n — 2) transverse coordinates, it is easy to see that the nontrivial components of the metric in this 
basis have the form 

77+- = ri-+ = 1, 1^'^ = -6'^ (5) 

so that the scalar product of two vectors can be written as 

A- B = A+B- + A-B+ - A- B (6) 

The momentum vector can also be written as — (p^ ,p), where p~ can be identified with the energy variable. 
In the light-front variables, the Einstein relation takes a linear form 

_ p^ +Tn^ 

This is a major difference from the conventional quantization on an equal-time surface. 

In the light-front variables, the Lagrangian density for a 0^ theory, for example, takes the form 

It is clear that the Euler-Lagrange equations, following from this, are only first order in the x^ derivative, which is a 
distinctive feature of light- front theories. The quantization of this theory has been discussed quite a lot in the literature 
and without going into details, we simply note here that the Feynman rules for this theory, at zero temperature, take 
the form 



2p^p — p^ — 




= -iXS''{p + q + r + s) (9) 
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where all the momenta, at the vertex, are assumed to be incoming. 

In going to finite temperature, as is well known, the interaction vertices of the theory are unaffected, but the 
propagators modify to reflect the periodicity (or anti-periodicity) of the field variables ^, ■ Let us generalize the 
theory in cq. (||) to finite temperature, following the conventional identification 

{0)fj = Tr e-l^"0 = Tr e-f^^' O (10) 

We can describe the resulting theory either in the real time formalism or in the imaginary time formalism and we 
discuss the two cases separately. 



A. Real time formalism 

Let us describe the propagators of the theory in the closed time path formalism |Q, ^ for simplicity. A similar 
structure for the propagators results in thermo field dynamics |l^], which we do not go into. It is well known that, 
in the real time formalism, the field degrees of freedom double and the propagators have a 2 x 2 matrix structure. 
In the case of the light-front scalar field theory with the conventional generalization in ([l^), the propagators for the 
doubled degrees of freedom have the forms 

'G++{p) = ^ + _ ' +27rnB{\p-\)S{2p+p- -LoD 

iG+_{p) = 2TT {e{-p-) + nB{\p-\)) 5{2p+p- 
iG^+{p) = 27r {e{p-) + nB{\p-\)) Si2p+p- ^LoD 

zG__(p) - - „ ^ _ ' , . +2nnB{\p-\)S{2p+p~ -cuD (11) 
2p+p — Lj^ — te ^ 

where we have introduced the bosonic distribution function 

1 



and have defined 



= p ^ + (12) 



It is worth remembering that Up involves only (n — 2) transverse components of the momenta, as opposed to the 
conventional theories where it depends on all the [n — 1) spatial components of the momentum. We also note that 
the ± subscripts refer to the "original" and the "doubled" degrees of freedom respectively. The propagators have the 
usual structure of a sum of the zero temperature part and the finite temperature part. However, the sign of trouble is 
already apparent in the form of the propagators in (^ij). The thermal distribution function does not seem to provide 
the necessary damping, as is usual in conventional thermal field theories. Namely, for p^ = Q (which is allowed by 
the delta function constraint), the distribution function diverges. As we will see soon, more difficulties arise in actual 
calculations. 

There are two kinds of vertices in the thermal field theory, "-I-" type and "— " type, with a relative sign difference 
between the two. However, at one loop, there is no mixing between the "original" and the "doubled" degrees of 
freedom. As a result, the one loop correction to the self-energy simply involves the tadpole graph (see Fig. 1), which 
can be readily evaluated. 

^^++(p) = -^/^^G++(fc) (13) 

Separating out the zero temperature part, the thermal correction to the self-energy can now be easily obtained. In 
fact, because of the delta function, the integral can be trivially done leading to 



iX 
2{2n) 



zHV^j(p) ^ / d^'knB{\k-\)5{2k-^k- -u^i) 



^d-^kT'^UBii^] (14) 



2(27r)"-i J Jo k+ \2k 




FIG. 1: One loop self-energy in (j)^ theory 

It is clear that the integral is divergent at the ultraviolet limit, fc+ oo, and needs to be regularized. Regularizing 
the power of k^ in the denominator yields 



mi^ip) = lim ^ / r-'k r ns ( ^\ 

++^^^ e^o 2(27r)"-i i io {k+Y+^ \2k+ J 



The form of the thermal correction in ( |15D is exact and is quite interesting. It shows that even though this repre- 
sents the thermal correction, it has a divergent part that is independent of temperature. Thus, it would seem that in 
a thermal background, the theory would require additional temperature independent counterterms beyond the ones 
needed for the regularization of the zero temperature theory. This is quite distinct from the behavior of the conven- 
tional thermal field theories and, if true, would cause enormous problems with the renormalizability properties of the 
light-front theories at finite temperature. Furthermore, since the integration over the (n — 2) transverse directions 
are yet to be performed, we see that the temperature dependent part of the amplitude is finite only in 1 + 1 dimen- 
sions. In any other dimension, the temperature dependent part diverges as well, requiring temperature dependent 
counterterms. We would like to emphasize here that, even though ( [l5| ) represents the temperature dependent part of 
the one loop self-energy, it does not vanish when temperature vanishes, namely as /3 — > oo. This is connected with 
the problem alluded to earlier, namely, since k~ can take a vanishing value, the limit /3 ^ oo of the distribution 
function is ambiguous (it is non-analytic at that point). Further problems arise when one studies the self-energy for 
the (j)^ theory on the light-front, but we will not go into the details of these here. These are serious problems which 
suggest that the naive generalization of ideas from conventional thermal field theories may not be appropriate in the 
case of light-front field theories. In the next section, we will analyze the source of these problems and propose the 
appropriate generalization for such theories. 



B. Imaginary time formalism 

In the imaginary time formalism, we rotate the theory to Euclidean space (imaginary time) and assume that the 
energy variable takes discrete values. Consequently, the energy integrals are replaced by a sum over the discrete 
Matsubara frequencies. In the light-front theories, this translates to replacing, in the scalar propagator, 

p~ — > 2innT (16) 

where T — ^ denotes temperature. As a result, the scalar propagator, in the imaginary time formalism, becomes 

^"^■^ 4i7rnrp+ — ojp ^ ^ 

The tadpole diagram. Fig. 1, is now straightforward to evaluate in the imaginary time formalism. 
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2(27r)"-i J ^ 4i7rnTfc+ 



n— — oo 
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Separating out the zero temperature contribution, we obtain, 



(18) 



4(27r)"-i y Jo V \^Tk+ J J 2(27r)"-i y Jo \'2k 

This is exactly the same expression as in (|l^) and, therefore, aU the subsequent analysis of the earlier subsection follows. 
The real time and the imaginary time formalisms give the same result which is, however, plagued by problems. We 
will discuss next the source of the problem. 

III. THE PROPER GENERALIZATION TO FINITE TEMPERATURE 

To understand the source of the problems in the last section, let us recapitulate briefly what happens in a conven- 
tional theory. In a conventional thermal field theory, the thermal part of the propagator represents the interactions 
of the particle with the thermal distribution of real particles in the medium. This is suppressed at high energies. In 
contrast, the propagators in eq. (|Tl|), as we have argued, do not provide the necessary damping. The form of the 
propagators are, of course, derived from the assumption that the ensemble average, in light-front theories, is given 
by eq. (|ic|). This is, in fact, where the problem lies. In a conventional theory, when one assumes that the ensemble 
average has the form 

{0)i3 = Tr e-/^" O, (20) 

it is understood that we are in a Lorentz frame where the heat bath is at rest. In fact, this is not a manifestly Lorentz 
covariant description. One can give a manifestly covariant description of thermal field theories |pl] , |l^ at the expense 
of introducing a velocity for the heat bath, u^, normalized to unity, namely 

u ■ u = u^u^ = 1 (21) 

and generalizing the ensemble average to 

{0)p = Tr e-^^-f" O = Tr e"'^""^" O (22) 
In a conventional thermal field theory, where the metric is diagonal and is of the form (+, —,—,•••,—), one can choose 



a rest frame of the heat bath corresponding to u^^ = (1,0,0, •__jO) consistent with (21) and, in this case, ( |22| ) would 
reduce to the conventional definition of ensemble average in (pO|). 

In contrast, a light-front description of a theory is manifestly relativistic. Intuitively, it is clear that it is not 
possible to have a heat bath at rest on the light-front. Eq. (10), on the other hand, is a generalization of the rest 
frame ensemble average ( pO| ) to light-front theories and, consequently, there are bound to be problems. Note that, for 
( p^ ) to reduce to (p^, we must have = (1, 0, 0, • • • , 0), which, with the light-front metric, gives 



2u u — u ■ u ~ 



This is inconsistent with (21) and this is another way of saying that we cannot have a heat bath at rest on the 
light-front. 

It is clear, therefore, that in dealing with light-front field theories, we must use a manifestly covariant description 
of the thermal field theories. With this, let us now discuss the proper generalization of real time and imaginary time 
formalisms for light-front field theories separately. In this section, we will restrict ourselves to scalar field theories 
and describe fermion theories in the next section. 

A. Real time formalism 

Let us assume that the heat bath is moving with a velocity u'^ subject to (|2l|). In that case, in the closed time path 
formalism, the propagators for the scalar field theory take the forms (for the doubled degrees of freedom) 

*G++(p) = -—-^^— + 27rnB{\u-p\)6{2p+p--Lul) 
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iG+-{p) = 2TT{0{~u-p) + nB{\u-p\))S{2p+p- -i^D 
iG^+ip) = 27r{9iu-p)+nBi\u-p\))S{2p+p- -L^^p) 

iG__{p) = , _ ' , - + 27rnB{\u-p\)S{2p+p- -L^D (23) 

Zp^p ~ LO^ — le ^ 

A simple choice for the velocity of the heat bath satisfying for example, is (in the light-front basis) 

u^ = -^(l,l,0,0,---,0) (24) 



in which case, we have 



u-p= ;^(P^ +P ) 



and the bosonic distribution function takes the form 

nB{\u-p\) = UB ^y^\P^ 

It is easy to see that this distribution function provides the necessary damping, on-shell, both at = and p+ — > oo. 

With this modification of the propagators, we can now re-evaluate the tadpole diagram (see Fig. 1). With (HJ), 
the temperature dependent part of the tadpole graph has the form 

lUi^lip) = -2(2^ / d-knB(^^\k+ + k-\y{2k+k--u;l) 

h'-^kT'^nBi^'^T^'] (25) 



2(27r)"-i J ^ Jo k+ ^\ 2V2k 

It is worth emphasizing that, unlike the corresponding expression with the naive generalization in ([l^), this integrand 
is well behaved in both the limits, fc+ and A:+ ^ oo as is expected of a thermal amplitude. As a result, it does not 
need any regularization. Furthermore, it vanishes at zero temperature, /3 —^ oo, as we would expect since it represents 
the thermal correction to the self-energy. However, in general, it cannot be evaluated in a closed form. In the high 
temperature limit, /3to ^ 1, and in four space-time dimensions, the integral has the value, 

^u'^_^lip)^-^ + OiPm) (26) 

There are several things to note from this result. First of all, this yields a thermal mass correction which, in the 
high temperature limit, has the form 

A^l = ^ = > (27) 

Namely, the thermal mass correction is positive as is the case in conventional theories. This is, in fact, crucial for 
restoration of symmetry at finite temperature. More interestingly, we note that the thermal mass correction coincides 
exactly with that obtained in a conventional thermal (scalar) field theory in four dimensions. 

Before closing this subsection on the real time formalism, let us note that the propagators ( |2^ ) satisfy the usual 
relations 

^G++(p)+^G__(p) - iG+-{p) + iG-+{p) 
iG++{p) - iG+-{p) = iG-+{p) - zG__(p) = iG^{p) 



{u • pY — • pY — a;^ + sgn(w • p) 

1 

'iG++{p) - iG-+{p) = iG+-{p)-iG — {p) = iG^ip) = 7 to — jz — to 5 ■ 1 \ (^8) 

[u ■ pY ~ [u ■ py — uj^ — le sgn(M • p) 

where Gr, Ga denote the retarded and the advanced propagators respectively and we have defined a vector which 
is orthogonal to (as well as to any vector in the transverse {n — 2) dimensional space) and has a space-like 
normalization, namely, 

u-u = {), u-u^-l (29) 
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For the choice in (H), uf" = ^(1, -1, 0, • • • , 0). 

For completeness, let us also note here the forms of the propagator in the formalism of thermo field dynamics in 
light-front scalar field theories, with proper generalization. 

iGnip) = „ , _ ' „ . +2TTnBi\u-p\)S{2p+p- ^ujD 
2p+p — u!^ + le ^ 

iGi2{p) = 2TmB{\u-p\)e^ 5{2p+p- -ujI) 

iG2i{p) = 2TinB{\u-p\)e^ 6{2p+p- -LoD 

zG22(p) = - + 27:71 b{\u ■ p\) 5{2p+ p- -uD (30) 

2p+p — uj^ — te 

Finally, let us note that since, in this case, we have two preferred vectors available, namely, and u^, any given 
vector can be uniquely decomposed as 

A^^iA- u)u'' - {A ■ uju" + At;, (31) 

where A!^ is transverse to both u, u. Similarly, any higher order tensor structure can also be decomposed with respect 
to these two vectors. In this way, a richer tensor structure arises in light-front theories at finite temperature than in 
conventional thermal field theories. 



B. Imaginary time formalism 

In the imaginary time formalism, in the covariant description, it is the variable {u ■ p) that is rotated to Euclidean 
space and takes discrete values. Thus, 

u ■ p 2i7mT (32) 

which, with the choice of (p^), leads to 

p~ ^ 2V2iTrnT - p+ = p~ - p^^ (33) 

where we have identified 

p- = 2V2iTmT (34) 

(An alternate way of doing the rotation is to decompose the momentum vector as in ( |3l| ) and rotate u-p while treating 
w - p as an independent variable.) As a result, the scalar propagator, in the imaginary time formalism, takes the form 

G(p) = —, ^rr— 7 = ^- (35) 

With this, let us calculate the tadpole diagram. Fig. 1, in the (f)^ theory. 



n(p) 



-— / d'-ifcV2r G(fc,n) 

' •' n=-oo 

27r)"-i J ^ AV2iTmTk+ - 



2(27r)«-i J ^^f-^ 4V2innTk+ ~ {loI + 2{k+Y) 

' fr-^k r^cothf-^+;(^"n (36) 



4(27r)"-i J Jo k+ \ AV2k+T 

Here, the factor of -\/2 arises from the Jacobian (because of the particular choice of the unit vector). Separating out 
the zero temperature part, then, leads to the thermal correction to the self-energy, 

_n.».M.^/.-*r^^„JMi:|^) (37, 



2(27r)«-i J Jo k+ \ 2V2k 



k+p 




P 



FIG. 2: One loop self-energy in cjy^ theory 

which coincides with the result calculated earlier in the real time formalism in (|2^) and all the subsequent analysis 
carries through. 

As another example, let us calculate the one loop scalar self-energy in a massive (jr" theory in 3 -|- 1 dimensions (see 
Fig. 2). In the imaginary time formalism, this has the form 

■,2 



nb) = 2(1^ / d^kV2TY,G{k,m)G{k+p,m) 

= / d3fc ^ 

2(27r)3i ^ 2{2^/2i^^nT - 



k+)k+ -Lol 2{2V2nrmT + p- - k+ - p+){k+ + p+) ~ ujI^j, 
The sum can be evaluated using standard formulae and leads to 
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{38) 



^^^^ 16(27r)3 J k+{k+ + p+) u.l+,+2{k++p+r _ c.g+2(fc+)^ _ 

2{k++p+) 2lT P 

X f coth + - coth i^'*-; ^"-^ ^ - ^1 ) (39) 

y 4\/2fc+r y 4V2(fc++p+)r 2V2TJJ 

If we use here the fact that p^ ~ 2\f2iT:nT as well as the periodicity of the hyperbolic function, the self-energy 
becomes (upon rotation to real time) 



(Pk 



— 2(fe++p+) 2FF yp +p ) 

( , L^2+2(fc+)2 , (ul +2{k+ +p+fW 

X coth '- ^^ — '— - coth — ^ (40) 

\^ 4x/2fc+T y 4V2(fc++p+)T yy 

It is now easy to take various limits of this expression. In fact, in the present case, we have more possibilities of 
taking limits than in a conventional thermal field theory [|2[ For example, we note that if we set p'^ = p = Q and 
take the limit p~ — s- 0, we obtain 

-n(p+ = 0,p- -> 0,p = 0) = (41) 
On the other hand, if we set p~ = p = and take the limit p+ 0, we obtain 

Finally, we can also set p~ = p+ = and take the limit p^O. In this limit, we obtain 

- n(p- . O.p- . 0,p ^ 0) . ^;r^l<fk ^c„.ed,^i4±|±l£ (43, 
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This shows that the three different ways of approaching the origin in the energy-momentum space lead to quite 
different results. Thus, light-front theories have a richer structure than the conventional thermal field theories also 
in this sense. Note, however, that as T ^ 0, all the three limits lead to a vanishing result, as would be expected in 
a zero temperature theory. Furthermore, an interesting question arises as to whether the three limits would lead to 
newer definitions of masses in light-front theories (in a conventional thermal field theory, we have only the screening 
mass and the plasmon mass corresponding to the two possible limits that are allowed). Even the question of what 
would correspond to the screening and the plasmon masses in such theories remains an open question. 



IV. FERMION THEORIES 



The fermion theories, on the light-front are more tricky simply because the number of degrees of freedom decreases 
in this case. Let us consider, for example, a free massive fermion theory on the light-front described by the Lagrangian 
density 

£ = i^ii-ff^d^ - to)?/' (44) 

We can, of course, add interactions, but, as we know, interaction vertices are not modified at finite temperature. 
Therefore, it is sufhcient to look at the free theory to determine the propagators at finite temperature. 
Let us define the light-front gamma matrices 

7± = -^(7"±7"-') (45) 

These are, in fact, nilpotent matrices, namely, 

(7^)' = (46) 

Defining the projection operators, 

^^^^ = ^7^7^-^(l±«n-i) (47) 
where a represents the Dirac matrices, it is easy to check that 

(PW)2=PW, p(+)p(-) =0 = P(-)P(+), PW+P(-' = 1 (48) 
With these projection operators, let us define 

^(±) ^ p(±)^ (49) 

Then, it follows from the properties of the gamma matrices that 

7"i/'W = = 7+V'(") (50) 
This allows us to write the Lagrangian density in the light-front variables as 



_i^(+)*^-(i-y . V + m)ip^~^ - ^V'^~-'*7'''(*7 • V + m)^'^"'"^ 



(51) 



It is clear now that only the 7/)'^+^* degrees of freedom are dynamical. The other degrees of freedom are related 
to these and can be eliminated. 

The fermion propagator, at zero temperature, has been derived long ago |^ and has the form 

^S,{x -v)^ (0|T+(V'(x)?(y))|0) = ^/2^ | ^e-^^<^-y^ {^^^^^ Jt) (^2) 

where "T+" denotes ordering with respect to Thus, we can identify 

MP) - V2 ( J^Z. -f,)-V-2 Z + r,. (53) 
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where p'^ denotes an on-shell momentum, namely, 



P , P = P, P = r, . (54) 



^f(p) = — T-— = -2—— (56) 



The second form of the propagator makes it very clear that, when properly normalized, it can be thought of as 
a projection operator and, consequently, its inverse does not exist. This is, therefore, not a suitable structure to 
generalize to finite temperature. On the other hand, the singular structure of the zero temperature propagator simply 
reflects the fact that there are constraints in the theory, namely, that not all the degrees of freedom are dynamical. If 
we eliminate the non-dynamical degrees of freedom, the entire theory can be recast in terms of i/j^^^ ■0'^+^^ variables. 
Therefore, the relevant propagator, from the point of view of the theory, is 

SF{x-y) - (O|r+(vW(x)0(+)'(y))|O) - (0|p(+)T+(V(x)V^(y))7"i^'+' |0) - J ^ e-^P^^-y^Spip) (55) 

This can, in fact, be calculated from ( ^ ) in a simple manner and is determined to be 

V2p+ _ V2p 
p'^ — rn^ + ie 2p^p^ — uj'p 

This is well behaved with the two point function, in this projected space, corresponding to 

p2 _ ^2 2p+p^ - UJp 
V2p+ ^ V2p+ 

and this propagator can be easily generalized to finite temperature. We note that this form of the propagator has 
also been calculated earlier directly from the field decomposition in p| . 

The fermion propagators at finite temperature, in the real time formalism (closed time path), now take the forms 

tS++{p) - V2p+( ' ~2TTnF{\u-p\)S{2p+p- ~u;l)) 

\ 2p+p — U!^ + te ^ ) 

iS+.ip) = -2V27rp+ {nF{\u ■ p\) ^ 0{~u ■ p)) 6{2p+p- -uD 

iS-+{p) = -2V2Trp+ {nF{\u-p\)^e{u-p))S{2p+p- -ljI) 



iS..{p) = V2p+{--— 2 --27TnF{\u-p\)S{2p+p- -u;l)) (57) 



where np represents the fermion distribution function 

nF{x) 



The propagators in thermo field dynamics can similarly be obtained and we do not go into this here. In the imaginary 
time formalism, the fermion propagator, for the independent degrees of freedom, takes the form 

S{p) = — = (58) 

2V2z(2n + l)^p+T - (oj2 + 2(p+)2) 



V. SUMMARY 



In this paper, we have described how light-front field theories can be generalized to finite temperature. We have 
shown that the naive generalization leads to problems and the origin of the difficulty is identified. Since light-front 
field theories describe relativistic systems, a covariant description of thermal field theories becomes necessary for 
the proper formulation of thermal light-front theories. We discuss scalar and fermion light-front field theories at 
finite temperature in detail, including issues such as non-analyticity of self-energy and tensor decomposition. Several 
open questions are also discussed. Light-front gauge theories at finite temperature as well as further applications are 
presently under study and will be reported later. 
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